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A numerical investigation of time-separated charge overlap measurements 
is carried out for the pion in the context of lattice QCD using smeared Wilson 
fermions. The evolution of the charge distribution function is examined and 
the expected asymptotic time behavior ~ g-C-Eg-m^)*^ where t represents the 
charge density relative time separation, is clearly visible in the Fourier trans- 
form. Values of the pion form factor are extracted using point-to-smeared cor- 
relation functions and are seen to be consistent with the expected monopole 
form from vector dominance. The implications of these results for hadron 
structure calculations is briefly discussed. 



The use of charge overlap distribution functions to examine quahtative 
features of lattice hadrons is now a well established technique [1-3]. How- 
ever, there is a need to go beyond qualitative measurements to investigate 
whether such methods are capable of extracting physical observables in the 
electromagnetic sector. For this reason the Fourier transform of the distribu- 
tion functions for tt and p mesons were examined in Refs. and 0. It was 
found that the equal-time functions contain a substantial contamination of 
intermediate states, preventing a reliable extraction of electric and magnetic 
form factors at available lattice momenta. 

Consider now the continuum Euclidean time-separated charge distribu- 
tion for zero-momentum pions: 
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J^''^ are the u, d flavor charge densities, and ' T' is time-ordering. Fig. 1 
shows a symbolic representation of the measurement in Eq. (|^). Assuming 
t > and inserting a complete set of states, this can be shown to result in 
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In the large Euclidean time limit, the sum reduces to a single term, and we 
find that 

QiUf,t) ^ ^-^fp^F^iq^) e-(^— (4) 
where (working in the flavor SU(2) limit) 

+ Jig 

Thus, by separating the currents in time it is in principle possible to damp 
out the intermediate state contributions and to measure form factors at arbi- 
trary lattice momenta. On a finite space-time lattice, however, achieving this 
elastic limit is problematical. The exponential damping factor, g-l^^-'^ir)*^ jg 
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not large and the fixed time positions of tlie initial and final pion interpolating 
fields limit the possible charge density separation times. Nevertheless, this 
paper demonstrates that the elastic limit can essentially be achieved at low 
momentum values on present-sized lattices. This has important implications 
for hadron structure function calculations on the lattice. 

We begin by examining the evolution of these distribution functions as the 
time separation, t, between the charge densities J^{x,t) and ^4(0) increases. 
On the lattice, we measure 

QiM h, ti, t) = E e^^'^Xt i2, ti, t), (6) 
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where we are now using lattice states and operators, and the fields 0*™' (t2) 
and 0^(0,0) are smeared and point pion interpolating fields, respectively. 
We use the exactly conserved lattice charge densities (which are non-local in 
time) for J^{x + r,ti), J^{x,t). Note that carrying out the sum on x in (|^) 
for all r takes a nontrivial amount of computer time, scaling as N^, where Ns 
is the number of spatial points in the lattice. One may show that for large 
time separations (^2 — ti, ti — t, t ^1) and in the continuum limit, 

QiMt,h,t2) ^ QtAf,h-t), Vt:{r,t,h,t2) ^ Pt:{r,t,-t). (8) 

The numerical results for this study are on 12 quenched j3 = 6.0 con- 
figurations of a 16^ X 24 lattice at k = .154 (the largest k value studied in 
Refs. and 0.) The pion interpolating sources were located at time slices 4 
and 21 and the time-separated current densities were positioned as symmet- 
rically possible in time between these sources. Two quark propagators per 
configuration, with origins at the positions of the interpolating sources, were 
necessary to reconstruct these relative overlap functions. Charge density 
self-contractions have been neglected in forming these quantities. 

Fig. 2 shows a log^Q plot of the spatial evolution of the charge overlap 
distribution function, 7^44, as the relative time separation between the d,u 
current densities is increased. It is seen that the distribution becomes flat- 
ter as the relative time is increased; there is also a substantial change in 
the shape. Points on the periodic lattice boundary are given by the filled-in 
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squares, which are seen to have their values raised. Note that these func- 
tions have been calculated using point-to-smeared correlation functions; the 
quarks are smeared over the entire volume of the lattice at time slice 21 using 
the lattice Coulomb gauged to produce zero- momentum quark and hadron 
fields. 

Fig. 3 represents a log^o'Plo^ of the Fourier transform of these distribution 
functions, Qf^, at the two lowest lattice spatial momentum values, |g| = ii/S, 
\/2(7r/8) function of relative time separation between current densities. 
The solid lines shown in this figure come from the expected asymptotic ex- 
ponential falloff specified by Eq. (^), using the (smeared) k = .154 data from 
Table 1 (m^r = .369, mp = .46) of Ref. [0, assuming the vector dominance 
form for the pion form factor: FT^{q^) = 1/(1 -|- q^/rn^p)- (We also assume 
the continuum relation Eq = (m^ + q'^Y^'^ .) Actually shown in this fig- 
ure are results for both point-to-smeared (O) as well as smeared-to-smeared 
(□) correlation functions. In all cases, the expected functional dependence 
~ Q-i^q-^^)* is present, indicating that by time step 7 or 8 single exponential 
behavior has emerged. This is similar to the number of time steps needed in 
hadron spectrum calculations. This behavior is remarkable because although 
we are damping out intermediate states as t increases, we are also moving 
closer to possible contaminations from the fixed interpolating fields at either 
time end. In fact, we do not see any indications of such contamination in 
the data. These correlation functions are also unusual because the asymp- 
totic approach is from below, indicating damping of negative terms in the 
(i, u correlation function. Ref. |^ indicates that these intermediate states, at 
least in the continuum limit, are primarily positive G-parity states^. 

Although the point-to-smeared and smeared-to-smeared results exhibit 
essentially overlapping error bars, the smeared-to-smeared values are sys- 
tematically low compared to the point-to-smeared results, indicating a slight 
dependence on the form of the interpolating field used. However, the results 
of Ref. indicate that such a dependence decreases as k — ^critical, that is, 
as the physical regime is approached. 

In looking back at Fig. 2, it is clear that the spatial distributions, Vf^{r), 
have substantial contaminations near the periodic lattice boundaries from the 
surrounding image sources . Why is it then that we have not attempted to do 
image corrections on the Fourier transform, Qf2{q), of these distributions? 
Consider an infinite spatial grid and associate the value of a function, F{x), 
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with each point x G {X} of the grid. Then, by doing a spatial translation 
on points outside the given primitive cell (of size L^) to equivalent points, 
associated with the same phase factor e*'''^, inside the cell x G {x}, it can be 
shown that 

^ e*«-^F(f) = J2 e''^'^F{x), (9) 

where for x G {x} 

= ^F(f+nL). (10) 

it 

{n = {nx,ny,nz), where n^^y^z are integers.) The terms with n ^ are the 
image contributions. Eq. (P) says that the discrete Fourier transform of F{x) 
on the infinite grid is the same as the transform of F{x) in the primitive cell 
for q values allowed in the original primitive cell. In our case this means 
there are no image corrections to Qf^ioJy assuming the distributions Vf^{f) 
are strictly periodic. 

By fitting time steps 7 through 10 of Fig. 3 with a single exponential of 
known slope and removing the kinematical factor in Eq. (^), a measurement 
of F.,^{q^) can be made. The results for the two calculated values using 
point-to-smeared correlation functions are shown in Fig. 4. (The error bars in 
Figs. 3 and 4 were determined from first and second-order single elimination 
jackknifes, respectively.) The values found are consistent within errors with 
vector dominance, shown as the solid line. Also shown in this figure are 
the results from a previous three-point-function simulation of the pion form 



factor] 10|. Comparison of the two types of results shows that the error bars 
of these different techniques are of the same order of magnitude for similar 
numbers of configurations. (Note the different (3 values and lattice sizes of 
the two simulations, however.) 

The significance of these results for hadron structure calculations goes 
beyond the specific findings discussed here for form factors. The ability to 
reach the elastic limit for hadron four-point-functions is crucial in attempt- 
ing to perform direct simulations of hadron structure functions. Up to the 
present, attention has been focused on the moments of such functions, which 
are given by the operator product expansion. These expansions are based 
upon separation of the short- distance physics, calculated perturbatively, from 
the long-distance part, evaluated on the lattice (or by other techniques) in 
the form of certain operator expectation values. The presence on the lattice, 
in the continuum limit, of all appropriate physical scales, combined with the 
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ability to make contact with the elastic limit of charge overlap four-point- 
functions, demonstrates that direct simulations of structure functions are 
feasible with current technology, at least at low lattice momenta. 

What additional quantities need to be calculated in order to initiate such 
a program? Clearly, the overlap function considered here, Qfl{q'^,t) is a 
sub-dominant quantity for structure functions; the short-distance physics 
will mainly be contained in the same- flavor Fourier transform: Q^iq'^it) 
or Q'^{q'^,t). Calculating this piece involves putting both electromagnetic 
current densities on the same quark line. This type of zero-momentum four- 
point function is difficult to simulate on the lattice, but can be obtained by 
combining fixed-source and quark-smearing techniques and is currently under 
investigation. In addition, the the current- current overlap distribution func- 
tion, (7r(0) I Ji(f, t)J;(0)|7r(0)), for both flavor-diagonal and nondiagonal 
currents is necessary in order to complete measurements of the two standard 
structure functions, Wi and W2- Of course it is the proton, not the pion, 
which is of most phenomenological interest. However, there is no barrier to 
applying these techniques to the proton as well, as long the lattices used are 
able to contain all the necessary length scales and the elastic limit in all these 
various four-point-functions can be demonstrated. 

To summarize, this paper shows that the elastic limit of rf, u charge over- 
lap distribution functions, in the context of the pion, can essentially be 
achieved at low momentum values on present-sized lattices. We do not claim 
that all the systematics of such measurements are yet completely under- 
stood, but the message that this limit is very close is extremely encouraging 
to attempts to provide direct calculations of structure functions from the 
fundamental QCD Lagrangian. 

This work was partially supported by the National Center for Supercom- 
puting Applications and utilized the NCSA CRAY 2 system at the University 
of Illinois at Urbana- Champaign. 
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Figure Captions 



Fig. 1: Symbolic representation of the d,u cliarge overlap measurement. 
The pion sources, shown as circles, are at fixed time positions in the 
lattice. 

Fig. 2: Evolution of the lattice d, u spatial charge density correlation func- 
tion, P44, for pions as the time separation, ti — t, between J|(f, ti), 
74(0, t) operators varies, as a function of non-equivalent lattice r val- 
ues, (a) ti-t = Q. {h) h-t = 3. (c) ti-t = Q. {d) ti-t = 10. 
Filled-in squares identify points on the periodic lattice boundary. 

Fig. 3: LogiQ plot of the Fourier transform of the spatial charge density cor- 
relation function, Q44, as a function of time separation ti —t between 
Jtifiiiji Jii^jt) operators. The upper points give the |g| = vr/S re- 
sults, the lower points correspond to |g| = \/2(7r/8). Results are given 
for both point-to-smeared correlation functions (O) as well as smeared- 
to-smeared correlation functions (□). 

Fig. 4: Pion form factor, F.^, as a function of g^/m^. The solid line is the 
monopole form from vector dominance. The data shown are the present 
results from point-to-smeared correlation functions (O) and previous 
results on the pion form factor from Ref. |0 (□). 
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